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摘 要 : 以 干旱 区 生态 脆弱 区 精 河 流域 为 研究 区 ,选取 1990 一 2016 年 气象 观测 日 值 数据 及 对 应 时 期 Landsat 系列 
影像 六 景 ,基于 地 表 能 量 平 衡 系统 模型 方法 ,对 1990 一 2016 年 蒸 散发 的 时 空格 局 进行 了 研究 。 结 果 表 明 :(1) 在 全 


球 升 温 的 大 背景 下 , 近 26 a 来 研究 区 燕 散发 整体 旦 增加 趋势 ,Morlet 小 波 分 析 显 示 存 在 5a.7 a、13 a 尺度 变化 周 
8], (2) 研究 区 燕 散 发 量 在 空间 呈现 南 高 北 低 特 征 , 藻 散发 减少 区 域 在 东北 未 利用 地 , 蒸 散 发 显著 增加 区 域 在 艾 比 
湖西 北部 和 南部 山区 。(3) 不 同 土地 利用 类 型 下 的 蒸 散发 量 由 高 至 低 为 :林地 > 水 体 > 草地 > 耕地 > 建设 用 地 > 未 利 


度 (T.)。 
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地 表 蒸 散发 是 陆 面 水 循环 中 最 重要 的 水 文 过 
程 之 一 ,也 是 联系 植物 气孔 行为 、 碳 交换 和 水 分 利 
用 的 关键 生态 过 程 ,同时 也 是 土壤 一 植被 一 大 气 联 
系 体 (SPAC, Soil-Plant- Atmosphere Continuum ) 中 水 
分 运动 的 重要 连接 桥梁 ""。 地 表 蒸 散 量 的 估算 不 仅 
是 流域 水 资源 评价 .管理 和 调控 的 重要 依据 ,而 且 
它 的 变化 还 影响 区 域 和 全 球 的 水 文 循环 ”。 

早 在 1802 年 Dalton 利用 包含 温度 .湿度 和 风速 
的 燕 发 公式 估算 地 表 蒸 散发 以 来 ,更 多 学 者 先后 提 
出 波 文 比 -能 量 平衡 法 ”、 蒸 散 空气 动力 学 方法 "、 
Penman-Monteith(P-M) 模 型 5 等 一 系列 蒸 散 发 观测 
法 和 计算 方法 。20 世 纪 70 年 代 后 ,伴随 着 遥感 技术 
的 快速 发 展 , 实 现 了 蒸 散发 研究 点 到 面 的 突破 ,并 
涌现 出 了 许多 燕 散 发 估算 的 空间 和 遥感 模型 ,如 地 表 
能 量 平 衡 指数 SEBI 模 型 "、 陆 面 能 量 平衡 算法 (SF- 
BAL, Surface Energy Balance Algorithm for Land ) "|、 
地 表 能 量 平衡 系统 (SEBS, Surface Energy Balance 
System) 等 。 其 中 ,SEBS 模 型 充分 考虑 干 湿 两 种 
极端 状况 ,尤其 是 针对 干旱 区 水 分 和 能 量 分 布 不 均 
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用 地 。(4) 与 地 表 藻 散发 表现 出 高 度 正 相 关 的 气象 因素 为 风速 ,地 表 参 数 为 高 程 (DEM) 与 地 表 比 辐射 率 (8); 与 地 
表 藻 散发 表现 出 高 度 负 相 关 的 气象 因素 为 近 地 表 空气 相对 湿度 ,地 表 参 数 为 温度 植被 干旱 指数 (TDVT) 与 地 表 温 


的 特点 ,能 在 具备 更 明确 物理 概念 的 同时 ,有 效 提 
高 蒸 散 量 遥 感 估 算 精 度 。 因 此 ,近年 来 在 国内 外 得 
以 较 广泛 应 用 ”。 

精 河流 域 作 为 艾 比 湖 流域 的 一 部 分 ,是 新 疆 生 
态 环境 最 脆弱 的 地 区 之 一 ,在 全 球 变 暖 的 大 背景 
下 ,多 年 径流 水 资源 不 断 减少 ,引发 一 系列 生态 环 
境 问题 ,精确 估算 区 域内 莹 散发 情况 可 为 优化 水 资 
源 配置 提供 依据 。 相 关 学 者 多 运用 参考 作物 计算 
法 研究 艾 比 湖 薰 散 量 的 时 空 变 化 特征 及 其 对 不 同 
气象 因子 的 敏感 性 ”"; 郭 家 新 等 "根据 水 热平衡 
原理 构建 TVETI 燕 散 指 数 反 演 艾 比 湖 流域 地 表 蒸 
散发 ; 张 晓 玉 等 3 结合 MODIS 数 据 和 气象 观测 数 
据 ,利用 SEBS 模型 估算 艾 比 湖 流 域 的 地 表 蒸 散 
发 。 但 对 其 子 流域 一 精 河流 域 的 蒸 散发 动态 变化 
及 其 机 理 研 究 较 少 ,本 文 借助 定量 遥感 手段 结合 
SEBS 模型 ,估算 精 河流 域 地 表 通 量 及 日 蒸 散发 量 ， 
运用 Morlet 小波 分 析 .重心 迁移 等 方法 分 析 蒸 散发 
的 动态 时 空 变 化 规律 和 特征 ,从 而 探索 该 模型 在 复 
林地 表 条 件 下 的 适用 性 。 
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1 研究 区 及 数据 


1.1 研究 区 概况 

Ti 19] jfi E (819 07' 52" -83*05' 48" E, 44?00'21" - 
45?00' 56" N ) hz FIRER SIKES PS TUR , E 
噶 尔 岔 地 西南 边缘 。 地 势 南 高 北 低 ,平均 海拔 320 
m, 夏 季 降 水 量 稀少 ,冬季 干燥 寒冷 ,大 陆 性 气候 特 
征 显著 。 区 内 风沙 天 气 多 ,持续 时 间 长 ,日照 充足 ， 
年 日 照 率 61% ~ 64% ,夏季 平均 日 照 时 数 在 10h 以 
上 ,干燥 的 气候 和 活跃 的 大 风 天 气 ,决定 了 研究 区 
内 强烈 的 蒸 散发 ,是 我 国 西北 干旱 区 生态 环境 严重 
退化 的 典型 区 域 之 一 (图 1)。 


Fig. 1 Location map of study area 
1.2 数据 来 源 
气象 资料 来 自 中 国 气象 局 气象 数据 中 心 网 站 
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(http://data.ema.en ) ,实际 观测 数据 选取 中 国 地 面 国 
际 交 换 站 精 河 气 象 站 (站 点 编号 :51334) 日 值 数据 ， 
时 间 起 止 为 1990 年 1 月 1 日 至 2016 年 12 月 31 日 。 

空间 遥感 数据 来 自 中 国 地 理 空间 数据 云 (http:/ 
www.gscloud.cn) ,分别 是 Landsat- 5 TM 的 1990 年 、 
1994 4. 1998 年 影像 数据 ;Landsat-7 ETM+ 的 2002 
年 .2007 年 .2011 年 影像 数据 ;Landsat-8 OLI. TIRS 
的 2013 年 .2016 年 影像 数据 ,月 份 均 为 9 月 (秋季 )， 
空间 分 辨 率 30 m, 轨道 号 :P146/R029, 云 量 均 低 于 
10%。 基 础 影像 预 处 理 及 空间 建 模 分 析 、 图 形 处 理 
分 别 用 ENVI 5.1、ArcGIS 10.1 软件 进行 操作 ,其 中 
基于 TM-destripe 工具 对 来 源 于 Landsat-7 ETM+ 的 
3 期 间 题 影像 进行 了 条 带 修复 处 理 。 


2 ” SEBS 模型 机 理 


基于 地 表 能 量 平衡 系统 (SEBS ) ,由 荷兰 学 者 Su 
于 2002 年 提出 外 。 模 型 理论 包括 三 个 部 分 :一 是 地 
表 物 理 参 数 反 演 ,主要 包括 地 表 反 照 率 、 地 表 比 辐 
射 率 、 地 表 温 度 、 归 一 化 植被 指数 (NDVI) .地 表 粗 烟 
度 ; 二 是 热传导 粗糙 度 长 度 运 算 , 主要 包括 参考 高 
度 下 的 气压 .温度 .湿度 和 风速 ,其 中 参考 高 度 用 大 
气 尺 度 气 象 模型 估算 ; 三 是 蒸发 比 计算 ,主要 包括 
向 下 的 太阳 辐射 和 长 波 辐射 (图 2)。 

模型 的 数学 基础 为 地 表 能 量 平衡 方程 号 
( 式 1): 


AE-R,-G-H (1) 
式 中 : AE 为 潜 热 通 量 (W.:m’) ;及 为 地 表 净 辐射 通 


图 2 SEBS 模 型 技术 流程 图 


Fig. 2 Technical flow chart of SEBS model 
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E(W-m”); H JJ dta E (W m? );C 为 土壤 热 通 量 
(W.m”)。 
地 表 净 辐射 通 量 的 计算 方法 "2( 式 2): 
R,=( - R, «(RI - R!) -(1 -eR (2) 
式 中 :及 为 地 表 净 辐射 通 量 (W.m2?) ;及 为 人 射 的 短 
波 辐射 量 (W.m”) ;a 为 地 表 反 照 率 ; R; 为 人 射 的 长 
波 辐射 量 (W.m”); R) 为 向 外 的 长 波 辐射 量 (W.m”); 
s 是 地 表 比 辐射 率 。 其 中 和 人 射 的 短波 辐射 是 利用 大 
阳 常 数 、 太 阳 入 射 角 、 相 对 日 地 距离 和 大 气 单 向 透 
过 率 计算 ;入 射 的 长 波 辐射 由 改进 的 史 蒂 夫 : 波 尔 
效 曼 定律 估算 而 来 ( 式 3); 向 外 的 长 波 辐射 采用 Ste- 
fan Boltzmann 方程 ,通过 地 表 比 辐射 率 和 地 表 温 度 
计算 ( 式 4)。 
Ri=1.08(-In7,) oT 


(3) 
R! - eg T. (4) 
式 中 : r, 为 大 气 单 向 透 过 率 ; o HBR WRZ 
REC [5.67x10* W * (m? KH ] ;7 为 空气 温度 (K); 
TT 为 地 表 温 度 (K)。 其 中 空间 空气 温度 数据 由 地 表 
温度 的 中 纬度 修正 函数 获得 。 
土壤 热 通 量 的 计算 方法 ( 式 5): 
G- RIF, *(1 -FF)T.-T) (5) 
式 中 : 工 .为 全 植被 覆盖 下 土壤 热 通 量 占 净 辐射 的 比 
(L.20.05) ; T, 裸 土 下 土壤 热 通 量 占 净 辐射 的 比 
( P, 20.315) ;了 为 植被 覆盖 度 , 通 过 归 一 化 植被 指 
数 (NDVI) 计 算 ( 式 6)。 
NDVI~NDYL (6) 
“~ NDVI,.. - NDVI,. 
感 热 通 量 计算 较为 复杂 ,是 把 地 表 参 数 和 大 气 
边界 层 参数 平均 风速 和 平均 温度 .地 表 位 温差 等 结 
合 , 通 过 整体 参数 化 求解 ,最 终 获 得 摩擦 风速 ws)、 


莫 宁 霍 夫 长 度 (2) , 感 热 通 量 (万 ) 。 
_H: z - d, e z - d, Zom 
p= eluf =) v [I Jd =] (1) 
3 
_ _ pe pð, 
L= ka (8) 


H - hy. pC, (0, - 0,) 


z-d, z-d, Zi "^ (9) 
PVC) 
式 中 :pw 为 摩擦 风速 (ms ) 让 为 卡门 常数 (=0.41); 


z 为 参数 测量 高 度 (m) ;qd 为 零 平 面 位 移 高 度 (m) 5 zo, 
为 动量 交换 粗糙 度 长 度 (m) ;zo 为 热量 交换 粗糙 度 


TEE (m) 记 为 莫 宁 - 奥 布 霍 夫 长 度 (m); 亚 ,和 亚 分 
别 为 动量 交换 和 热量 交换 MOS 稳 定 度 订正 的 数 ;p 为 
空气 密度 (kg"m”);0, 为 空气 定夺 比热容 (J'kg K) 
(C,21005) ;9。 和 9, 分 别 为 地 表 位 温 (K) 和 参数 测量 
高 度 z 处 的 大 气 位 温 (K)。 
蒸 散 比 的 计算 过 程 中 假设 : 
CD 假设 在 极端 干旱 状态 下 ,地 表 无 蒸 散 发 , 湾 
热 通 量 为 0, 则 : 
H,-R,-G (10) 
@ 假设 在 极端 湿润 状态 下 ,水 分 充足 , 感 热 通 
量 趋 近 于 0, 此 时 蒜 散 发 速率 最 大 , 则 : 
H,-R,-G-AE, (11) 
e pu: 
"^ kg0.61(R,-G)AÀ 


_ 1 1z-do1 z-d, E 
x. P ) v( L jv) (13) 


HED O, RHE FEARI RO 
(3X 15): 


(12) 


(14) 


- AE .AAE, 
R,-G R,-G 
式 中 : A ,为 相对 蔡 散 比 ;也 为 极端 干旱 条 件 下 的 感 
热 通 量 ;及 ,为 极端 湿润 条 件 下 的 感 热 通 量 ; AE, 为 
湿润 条 件 下 的 潜 热 通 量 ;L, 为 极端 湿润 条 件 下 奥 布 
霍 夫 长 度 ;r, 为 极端 湿润 条 件 下 外 部 阻抗 。 

最 终日 燕 散 发 估算 借助 蒸 散 比 完成 : 


Ro 
ET, =8.67x 10 x A, x —"- 
， x 


IP : ETa H HRE; An 为 瞬时 蒸发 比 , An = 
A ;Ra 为 日 净 辐 射 通 量 ; A 为 气 化 潜 热 ; p 为 大 气 

SEBS 所 采用 的 动量 粗糙 度 模 型 适用 于 分 布 均 
匀 的 低 矮 植 被 。 因 地 区 土地 覆 被 类 型 等 条 件 的 不 
同 ,会 导致 地 表 粗 糙 度 存在 明显 差异 ,进而 对 模型 
模拟 精度 产生 较 大 影响 。 本 研究 根据 研究 区 植被 
覆盖 的 实际 情况 ,基于 二 阶 闭合 理论 ,对 动量 粗糙 
BEBE zon ) 零 平面 位 移 ( d, ) 的 计算 进行 调整 ,以 
减少 因 地 表 植被 分 布 不 均 引 起 的 蒸 散发 反 演 的 不 
确定 性 。 其 中 动量 粗糙 度 长 度 的 计算 区 分 部 分 植 
被 覆盖 和 完全 植被 覆盖 两 种 情况 ( 式 17); 零 平面 位 
移 借助 叶 面积 指数 计算 ( 式 18)。 


(15) 


(16) 


202103.00099v1 


chinaXiv 


za * 0.3h(C, LAI)", (0 LAI «I) 


i 17 
oaii : %) (1 <LAI<7.5) M 
d,=1.1 Infi + (CLAD | (18) 


式 中 :hh 为 冠 层 高 度 ; zo 底层 粗糙 长 度 , 取 0.01; C, 
为 叶片 拖 搜 系数 , 取 0.2;LAI 为 叶 面积 指数 ,可 利用 
植被 覆盖 度 获取 。 
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3 ”精度 验证 


对 于 本 研究 模型 结果 的 验证 ,主要 采用 “点 ”和 
“ 面 ””“ 相 结合 的 方法 。“ 点 "验证 主要 通过 精 河 
气象 站 点 的 日 实测 蒸发 数据 完成 ( 表 1)。 其 中 ,对 
比 蒸发 血 数 据 ,SEBS 模型 菊 散 量 的 平均 误差 率 
为 -8.39% ; 对 模型 中 地 表 参 数 的 验证 ,地表 温 度 .大 
气温 度 的 平均 误差 率 分 别 为 11.22% 12.79%. 


表 1 SEBS 模型 精度 验证 结果 


Tab.1 SEBS model accuracy verification results 


日 期 Ea E SEBS 相对 实测 地 表 遥感 反 演 相对 实测 大 气 MARS TR 相对 

(E601B)/mm ” 蔡 散 量 /mm ”误差 /% 温度 /%C 地 表 温 度 /"C ”误差 /% 温度 /人 C 大 气温 度 /CC 误差 /% 
1990-09-05 2 1.706 -14.72 17.3 19.939 13.23 15.6 14.647 -6.11 
1994-09-30 24 2.099 -0.03 18.5 19.822 6.67 6.2 13.74 -15.19 
1998-09-25 2.8 2.336 -16.56 19.8 22.771 13.05 18.2 14.621 -19.66 
2002-09-28 2.8 2.741 -2.13 19.1 23.265 17.9 17.6 16.815 -4.46 
2007-09-10 4.1 3.701 -9.73 2731 25.254 -9.68 21.5 20.235 -5.88 
2011-09-05 4.4 4.254 -3.31 28.1 24.329 -15.5 23.5 19.432 -17.31 
2013-09-18 3.5 3.48 -0.59 18.5 19.776 6.45 16.6 22.158 25.08 
2016-09-25 2.7 2.389 -11.52 29.12 31.406 7.28 23.65 25.689 8.62 


“ 面 " 的 验证 过 程 中 ,由 于 Penman-Monteith 模 
型 具有 很 好 的 物理 基础 ,作为 统一 的 标准 计算 方 
法 ,无 需 进 行 地 区 校正 和 改变 任何 参数 即 可 适用 于 
世界 各 地 ,上 且 具 有 和 较 高 精度 和 良好 可 比 性 “” , 故 在 
参考 “点 ”验证 结果 的 基础 上 ,将 Penman-Monteith 
模型 SEBS 模 型 反 演 的 区 域 整体 蒸 散发 结果 进行 对 
比 验 证 ,结果 表明 SEBS 模 型 呈现 的 多 年 趋势 与 Pen- 
man- Monteith 模型 数据 较为 一 致 (图 3) ,能 够 很 好 
的 适用 于 精 河 地 区 。 


4.5r —#— P-M —— SEBS 


40b o 线性 (P-M) ----- 线性 (SEBS) 


LO Wer uA y=0.2277x+0.4247 
z R?-0.5871 


0 n n n n n 1 " 2 
1990 1994 1998 2002 2007 2011 2013 2016 
年 份 


图 3 模型 精度 验证 
Fig. 3 Accuracy verification of model 


4 KARRETE 


4.4 蒸 散发 时 间 变 化 
研究 区 模型 反 演 日 蒸 散 量 在 26 a 间 整体 呈 增 加 


趋势 ,表现 为 1.035 mm- (103) ';. Ht} , 1990 年 至 
1998 年 变化 浮动 小 ,均值 分 别 为 1.5 mm 1.647 mm, 
1.41 mm; 在 1998 一 2011 年 有 显著 增加 ,均值 增加 
1.89 mm; 主 要 体现 在 艾 比 湖 湖区 .南部 及 东南 部 山 
区 、 精 河流 域 中 上 游 平 原 区 。2011 一 2013 年 减少 ， 
均值 减少 0.821 mm;2013 一 2016 年 再 次 增加 ,均值 
增加 1.316 mm( 图 4)。 

为 了 解 研究 区 蔡 散 发 变化 的 周期 特征 ,利用 精 
河 气象 站 (站 点 编号 :51334)1990 一 2016 年 历史 观 
测 数据 ,计算 Morlet 小 波 系 数 并 绘制 小 波 方 差 图 及 
小 波 系数 实 部 等 值 线 图 。 小 波 方差 图 显示 ,在 
1990 一 2016 年 研究 期 内 有 3 个 峰值 ,分 别 在 5 a、7 a, 
13 a 的 时 间 尺 度 上 (图 5)。 其 中 ,13 a 大 尺度 周期 对 
照 图 6 周期 震荡 特征 非常 明显 ,经历 了 7 个 正 负 循 
环 交 蔡 过 程 ,日 2014 年 出 现 的 正 等 值 线 至 2016 年 
尚未 闭合 ,表明 未 来 一 段 时 间 内 蒙 散 量 将 会 偏 高 ， 
既 印 证 了 前 述 研 究 期 内 的 蒸 散发 特征 ,又 表现 出 研 
究 区 未 来 蒸 散发 持续 增高 的 趋势 。7 a 时间 尺 度 上 
的 小 波 方差 值 略 高 于 5a 时 间 尺 度 , 从 图 6 中 可 以 看 
出 ,7a 尺 度 的 周期 变化 下 同时 般 套 了 5 a 的 周期 震 
i ,说 明 精 河流 域 蒸 散 发 存在 不 同 尺度 周期 的 复杂 
髓 套 结构 , 即 小 尺度 变化 和 藤 套 在 大 尺度 的 变化 下 ， 
这 与 其 他 学 者 研究 一 致 2。 
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(a) 1990-09-05 (b) 1994-09-30 


(d) 2002-09-28 (e) 2007-09-10 


(g) 2013-09-18 (h) 2016-09-25 


(c) 1998-09-25 


(f) 2011-09-05 
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Fig.4 Spatial distribution of evapotranspiration in remote sensing inversion in the study area 


小 波 方差 


0 2 4 6 8 10 12 14 16 
时 间 尺 度 /a 
图 5 小 波 方差 


Fig.5 Wavelet variance diagram 


42 土地 利用 类 型 影响 下 的 蒸 散 发 空间 特征 

研究 区 藻 散 发 空间 分 布 不 均衡 ,湖泊 河流 上 
下 游 \ 山 地 的 蒸 散 量 差 异 明显 ,基本 呈现 南 高 北 低 
的 空间 分 布 特征 , 即 南 部 林地 出 现 相对 高 值 ,北部 


1 —. 
1990 1995 2000 2005 2010 2015 
年 份 
图 6 小 波 系数 实 部 等 值 线 


Fig. 6 Real coefficient contour of wavelet coefficients 


除 艾 比 湖区 为 高 值 区 外 , 干 裸 的 未 利用 地 出 现 最 低 
值 。 由 图 7 可 知 , 燕 散 发 增加 及 基本 不 变 区 域 所 占 
面积 最 多 , 达 74.54% , 主要 集中 在 北部 湖区 .中 部 绿 
洲 区 及 南部 山区 ; 蒸 散 发 减少 区 域 主要 集中 在 东北 
部 未 利用 地 ,所 占 面积 为 25.46%。 
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图 7 近 26a 研 究 区 蒸 散 发 空间 变化 趋势 
Fig. 7 Trend of evapotranspiration in the near 26 years study 


area 


参考 中 国土 地 资源 分 类 系统 六 ,对 研究 区 进行 
土地 利用 分 类 。 在 借助 重心 迁移 模型 明确 了 人 研 
究 区 近 26 a 不 同 土地 利用 类 型 的 空间 重心 迁移 特征 
基础 上 (图 8) ,对 比分 析 了 研究 区 莹 散发 空间 变化 
趋势 (图 7)。 结 果 表 明 : 近 26 a, 耕 地 、 林 地 .草地 均 
呈现 向 东 的 迁移 趋势 ,耕地 整体 向 东 迁 移 4.99 km, 
而 林 .草地 呈现 向 东南 方向 的 迁移 趋势 ,分 别 迁 移 
13.96 km 、10.37 km ,使 蒸 散 发 轻微 增加 区 域 集中 分 
布 在 艾 比 湖 湖 区 及 其 以 南 的 水 域 充沛 区 域 。 蒸 散 
发 呈现 高 值 的 水 体 ,向 西北 方向 迁移 7.21 km, 因 而 
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图 例 使 蒸 散 发 显著 增加 的 区 域 主 要 出 现在 艾 比 湖西 北 


部 。 蒸 散发 呈现 低 值 的 建设 用 地 .未 利用 地 有 明显 
的 分 异 。 其 中 ,由 于 绿洲 城市 的 扩张 ,建设 用 地 整 
体 向 西南 方向 迁移 23.67 km ,未 利用 地 呈现 向 东北 
方向 的 迁移 趋势 , 共 迁 移 3.08 km, 导 致 莹 散发 减少 
区 域 主要 发 生 在 东北 未 利用 地 和 中 部 托 里 乡 旦 达 


荐 沙漠 。 
5 驱动 因子 分 析 


5.4 自然 驱动 因子 

蒸 散 发 作为 表征 气候 的 因子 之 一 ,其 变化 决定 
于 与 其 相关 的 气象 因子 的 共同 作用 结果 。 为 定量 
分 析 各 个 气象 因子 对 燕 散 发 变化 的 贡献 ,对 气温 、 
相对 湿度 、 降 水 量 、 风 速 、 太 阳 辐 射 、 日 照 时 数 及 日 
较 差 进 行 了 相关 性 分 析 ( 表 2)。 

针对 研究 区 地 表 薰 散发 而 言 , 由 于 研究 区 西北 
部 阿拉 山口 的 存在 ,风速 与 其 有 着 最 高 的 正 相 关 
性 ,而 近 地 面 空气 相对 湿度 则 对 其 表现 出 最 高 的 负 
相关 影响 。 需 要 强调 的 是 ,气温 相 较 于 其 他 气象 因 
子 与 研究 区 蒸 散发 正 相 关 性 的 显著 水 平 虽然 较 低 ， 
但 由 于 其 决定 着 空气 中 饱和 水 汽 含量 和 水 汽 扩 散 
的 速度 ,因此 对 蒸 散发 的 直接 影响 仍 不 可 忽视 。 联 
合 国政 府 间 气 候 变 化 专门 委员 会 (IPCC) 的 报告 显 
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图 8 近 26 a 人 研究 区 各 土地 利用 类 型 空间 重心 迁移 


Fig. 8 Migration of the spatial center of gravity of various land use types in the near 26 years study area 
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表 2 实际 蒸 散 量 与 气象 参数 相关 性 分 析 


Tab.2 Correlation analysis between actual evapotranspiration and meteorological parameters 


相对 湿度 降水 量 


气温 


风速 太阳 辐射 日 照 时 数 


日 较 差 


实际 蒸 散 量 0.634" -0.796* -0.337" 


注 :*、** 分 别 表示 在 0.05 和 0.01 水 平 上 显著 相关 。 


zR 7" ,2015—2019 年 ,全 球 平均 气温 较 工 业 化 前 时 
代 升 高 了 1.1 C ;对 研究 区 而 言 , 近 60 a 的 平均 气温 
以 0.33 % (10a 六 的 气候 倾向 率 呈 显著 (ao=0.05) 上 
升 趋势 ”…。 而 反 演 结果 显示 出 蒸 散 量 在 近 26a 的 持 
续 增 加 ,在 一 定 程 度 上 响应 了 气候 变 暖 的 事实 。 在 
2013 年 出 现 的 低 值 , 则 主要 是 由 于 影像 日 期 前 1 d 
的 天 气 过 程 导致 气温 与 风速 均 较 低 所 致 。 
Ut dei BE CT.) He EE ET SS Ce) .地面 数 
字 高 程 (DEM) 以 及 温度 植被 干旱 指数 (TVDI)4 个 地 
表 参 数 ,与 研究 区 蒸 散发 进行 相关 性 分 析 ,结果 表 
HH .DEM、 地 表 比 辐射 率 与 藻 散 发 呈正 相关 关系 。 
其 中 海拔 超过 500 m 后 , 除 在 2800 m 人 处 由 于 天 山 云 
杉 分 布 海拔 界线 会 出 现 小 幅度 下 降 外 , 蔡 散 发 整体 
随 DEM 增 加 呈 显 著 上 升 。 加 之 ,地 形 能 够 影响 地 表 
辐射 能 量 的 吸收 状况 , 故 蒸 散 量 也 随地 表 比 辐射 率 
的 增加 而 增加 。 由 于 研究 区 未 利用 地 多 为 岩石 X 
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BE .盐碱地 ,植被 覆盖 率 少 ,地 表 升 温 快 ,因此 地 表 
温度 与 蒸 散发 旺 负 相 关 关 系 。 而 TVDI 与 蒸 散发 的 
负 相 关 性 ,表明 了 表层 土壤 水 分 含水 量 对 莹 散发 的 显 
著 影 响 , 即 士 壤 水 分 含水 量 越 少 , 蔡 散 发 越 小 (图 9)。 
5.0. 人 为 影响 因子 

不 同 的 土地 利用 会 使 土壤 湿度 和 地 表 温 度 状 
况 发 生 改变 ,进而 对 区 域 蒸 散 发 产生 影响 。 由 表 3 
可 见 , 近 26a 研 究 区 建设 用 地 水体、 耕地 整体 呈 增 
加 趋势 ;林地 .草地 .未 利用 地 整体 呈 减 少 趋势 。 其 
中 ,在 人 口 不 断 增 长 及 绿洲 农业 的 发 展 影响 下 , 耕 
地 面积 增长 迅速 ,动态 度 为 12.24% ; CHE Wb R 5 c 
少 趋 势 ,但 自 研 究 区 实施 “三 北 ” 防 护林 、 退 耕 还 林 、 
等 林业 重点 工程 及 草原 生态 保护 奖 补 政策 以 来 , 面 
积 逐 渐 恢 复 。 因 此 ,致使 林地 .草地 ,耕地 成 为 蒸 散 
发 高 值 区 。 由 图 10 可 以 看 出 ,牲畜 数量 随 人 口 总 数 
增加 ,也 呈现 增加 趋势 。 牲 冀 数 量 作为 衡量 草场 合 
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Fig. 9 Scatter plot of evapotranspiration and surface parameters 
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表 3 1990 一 2016 年 研究 区 土地 利用 动态 度 变化 


Tab.3 Dynamic changes of land use in the study area from 1990 to 2016 [96 
年 份 林地 水 体 草地 耕地 建设 用 地 未 利用 地 
1990—1994 -0.50 1.34 10.32 2.17 13.48 -1.62 
1994—1998 -7.714 0.37 -9.10 14.58 1.71 2.07 
1998—2002 -7.59 10.95 12.43 12:11 1.40 =2.59 
2002—2007 —0.88 -8.18 -1.54 7.24 16.21 0.66 
2007—2011 1.24 1.50 2.82 1.06 6.61 -0.38 
2011 一 2013 27.49 -5.26 7.18 10.40 9.49 -2.83 
2013—2016 36.85 15.31 -6.43 -L51 14.36 1.62 
1990—2016 -0.82 0.99 -0.22 12.24 18.89 -0.33 
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图 10 研究 区 总 人 口 与 牲畜 数量 的 关系 


Fig. 10 The relationship between the total population 


and the number of livestock in the study area 


理 利用 的 重要 指标 ,其 增长 的 快慢 可 直接 导致 草场 
面积 变化 ,说 明 随 着 人 类 对 大 自然 改造 能 力 的 大 幅 
度 提 高 ,能 够 直接 影响 下 热 面 土地 歼 被 的 变化 , 进 
而 在 一 定 程 度 上 影响 下 垫 面 的 蒸 散 发 变化 。 


区 周边 植被 稀少 、 未 利用 地 面积 增加 ,使 得 研究 区 
蒸 散 发 减少 区 域 出 现在 东北 未 利用 地 。 但 因 2007 
年 艾 比 湖 国 家 级 自然 保护 区 的 成 立 及 研究 区 统一 
管理 .优化 水 资源 ,实施 节 水 灌溉 工程 , 自 2016 年 
起 , 艾 比 湖 湖面 开始 恢复 ,至 2017 年 有 328 km 湖面 
面积 扩充 ,由 此 使 得 研究 区 显著 增加 的 蒸 散发 区 
域 出 现在 了 保护 区 西北 部 。 本 文 研究 结论 如 下 : 
(1) 在 全 球 气候 变 暖 的 背景 下 , 精 河流 域 
1990 一 2016 年 蒸 散发 整体 呈 波 动 性 增加 趋势 ,随时 
间 推 移 存在 着 Sa 7 a 13 a 的 三 类 尺度 的 周期 变化 


规律 。 在 空间 上 基本 呈现 南 高 北 低 的 特征 ,26 a 间 
6 结语 蒸 散 发 减少 区 域 在 东北 未 利用 地 , 蒸 散 发 显著 增加 


区 域 获 散发 受 自然 和 人 为 因素 共同 作用 的 影 
响 ,在 气候 变化 导致 全 球 升温 的 背景 下 ,我 国 西北 
诸 河 流域 整体 蒸 散发 增 大 ,本 文 研究 的 结果 同样 
显示 了 蒜 散 发 显著 增加 的 事实 ,并 且 其 时 空 分 布 格 
局 也 在 发 生 着 变化 。 随 着 区 域 经 济 的 发 展 和 城市 
化 进程 加 快 ,虽然 越 来 越 多 的 自然 生态 下 垫 面 被 人 
工 下 垫 面 所 取代 ,天然 渗 透 层 被 不 透水 层 替 代 ,但 
对 营 散 发 贡献 较 少 的 城市 用 地 仍 不 可 忽视 。2012 
年 以 来 ,研究 区 在 践 行 生 态 文明 和 绿色 理念 的 过 程 
中 新 增城 市 绿地 面积 0.127 km ,使 得 建设 用 地 蒸 散 
发 逐年 增加 。 

此 外 还 需 关 注 的 是 ,虽然 26 a 间 研究 区 蒸 散发 
整体 呈正 向 发 展 ,但 在 空间 上 不 均衡 ,这 种 变化 在 


区 域 在 艾 比 湖西 北部 和 南部 山区 。 

(2) 在 研究 区 域内 ,不 同 土地 利用 类 型 下 的 落 
散发 表现 为 :林地 > 水 体 > 草地 > 耕地 > 建设 用 地 > 未 
利用 地 ,说 明 区 域 获 散 量 的 大 小 不 仅 受 区 域 气候 、 
自然 地 理 等 条 件 影响 ,还 与 下 热 面 土地 利用 、 植 被 
盖 度 等 有 关 。 随 着 经 济 的 发 展 ,自然 生态 下 热 面 被 
人 工 下 执 面 所 取代 的 同时 ,也 在 影响 并 改变 着 区 域 

(3) 研究 区 风速 与 地 表 蒸 散发 具有 显著 正 相 关 
性 , 近 地 面 空气 相对 湿度 与 地 表 蒸 散发 具有 显著 负 
相关 性 ;地表 参 数 相关 分 析 表 明 ,地 表 比 辐射 率 
(e) DEM 与 蒸 散 发 具有 显著 的 正 相 关 性 ,地 表 温 度 
(T) 温度 植被 干旱 指数 (TVDDU) 与 蒸 散发 具有 显著 
负 相 关 性 。 
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Evapotranspiration in the Jinghe River Basin based on the 
surface energy balance system 
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Abstract: The arid ecological fragile area of Jinghe Basin was selected as the research area. Daily meteorological 
observation data from 1990 to 2016 and six scenes of corresponding Landsat series images were investigated. The 
results show that under the context of global warming, evapotranspiration in the study area as a whole has increased 
over the past 26 years and Morlet wavelet analysis shows that there are 5 a, 7 a, 13 a scale change cycles. 
Additionally, evapotranspiration in the study area is spatially high in the south and low in the north. The 
evapotranspiration was reduced in the unused land in the northeast and significantly increased in the northwestern 
part of Lake Aibi and the southern mountainous area. Evapotranspiration under different land-use types followed the 
pattern from high to low of forest land > water body > grassland > cultivated land > construction land > unused land. 
Among the measured meteorological factors, wind speed was highly positively correlated with surface 
evapotranspiration, along with elevation and specific surface emissivity. The meteorological factors that show a 
highly negative correlation with surface evapotranspiration were near-surface air relative humidity, and the surface 
parameters were temperature vegetation drought index and surface temperature. 


Keywords: SEBS model; evapotranspiration; space-time pattern; Jinghe River Basin; Xinjiang 


